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ABSTRACT 


Vertical velocities have been computed for the lower stratosphere for two independent winter periods, by em- 
ploying a form of the adiabatic method. The regions studied were in both cases outside the polar vortex. The 
flow pattern was divided into stationary (long-wave) and transient (short-wave) components. The vertical 
velocity pattern associated with the stationary long wave is precisely that described by Kochanski [3] ; ie., the air 
rises in moving from warm troughs to cold ridges. The pattern associated with the short waves is more com- 
plicated. There is a maximum of warm advection in the vicinity of short-wave ridges, and cold advection near 
troughs. Local temperature changes, however, very nearly compensate the advection, with the net result that in 
the mean the vertical velocities associated with short-wave patterns are small, but tend to be positive near ridges 
and negative near troughs. Superimposition of the short and long wave, however, can lead to any conceivable com- 
bination of signs of advection, local temperature change, vertical velocity, and position with respect to ridge or 


trough. 


The single parameter which is most useful in specifying the vertical velocity is the temperature advection. 


1. INTRODUCTION 


In the middle 1940's, spurred by recent developments 
in atmospheric sounding techniques, meteorologists began 
to investigate in considerable detail the nature of the 
vertical motions in the lower and middle troposphere and 
their relation to the horizontal flow. Now, in the late 
1950's, as further improvements have increased the den- 
sity and reliability of data in the lower stratosphere, it is 
feasible to extend our investigations to those levels. The 
present study is offered as an early step in that direction. 
A hypothesis will be formulated concerning the manner 
in Which the vertical velocities are related to the troughs 
and ridges in the lower stratosphere, in the hope that this 
may serve as a basis for further studies. 

The earliest attempt at an empirical evaluation of the 
vertical velocities in the lower stratosphere was that of 
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Kochanski in 1951 (cf. [3]). Mean soundings at troughs 
and ridges were compared, and vertical velocities were 
computed on the assumption of adiabatic motion by com- 
puting the time it would take for parcels of air to move 
along a contour from ridge to trough. Average vertical 
velocities of the order of 3 cm. sec." were found near the 
tropopause, diminishing to almost negligible amounts 
just above 100 mb. Brockmeyer [1] studied vertical 
velocities in the lower stratosphere for the period Janu- 
ary 1-10, 1956 by means of the usual adiabatic method. 
His results were quite similar to those of Kochanski; 
i.e., values of vertical velocity ranging up to 3 to 4 cm. 
sec. at 200 mb. with smaller values at higher levels, and, 
in general, subsidence from ridge to trough and rising 
motion from trough to ridge. 

Craig and Hering [2], in a study of the sudden warm- 
ing in the polar stratosphere in January 1957, estimated 
that to account for changes in temperature along a tra- 
jectory, assuming adiabatic processes, vertical velocities 
of about 5 em. sec.’ were required at 25 mb. This, of 
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course, Was a somewhat anomalous situation. Teweles 
[4], in discussing the same situation, estimated average 
rates of subsidence of 1 to 3 cm. sec.* over somewhat 
longer paths. 


2. METHOD 


The technique chosen for computing vertical velocities 
is a modified version of the adiabatic method. This 
method is particularly well suited to the stratosphere 
because of the large stability in that region. The basis 
of the adiabatic method is the equation: 


(1) 


The assumption is made that the horizontal advection of 
temperature (V-V7') is very nearly equal to constant pres- 
sure advection. 

To simplify the computations it is also assumed that 
the lapse rate (y) may be neglected in the stratosphere, 
in comparison with the adiabatic lapse rate (yaa). This 
is not likely to be in error, in the regions studied, by more 
than 15 or 20 percent. Finally, then, if the vertical 
velocity (w) is measured in units of cm. sec., and if the 
units of advection and local temperature change are 
°K.(12 hr.)~, equation (1) may be rewritten 


(2) 


In most previous applications the right hand side of 
equation (2) has been evaluated by estimating tempera- 
ture changes along trajectories. This is unsatisfactory 
for two reasons. First, there is the difficulty of construct- 
ing trajectories from synoptic maps; second, the vertical 
velocities computed are then averages both in time and 
space over the length of the trajectory. Since a parcel of 
air may travel hundreds or even thousands of kilometers 
in 12 hours, this leads to considerable ambiguity in inter- 
preting the results with respect to the synoptic pattern. 

By evaluating the two terms on the right hand side of 
equation (2) separately, a method is developed by which 
adiabatic vertical velocities, although averages over 12- 
hour time intervals, refer to specific areas, which, in the 
present case, are 10° latitude by 10° longitude. There 
is no longer any question as to where, along a trajectory, 
the vertical velocities occurred. 

The basic grid is indicated in figure 1. The input in- 
formation was the temperature at every point in the 
5° x 5° grid (numbers in fig. 1) and mean winds within 
every 5° box (indicated by letters). The temperatures 
were estimated to the nearest 14° K. on the basis of an iso- 
thermal analysis made for every 2° K. The vector winds 
were estimated subjectively to the nearest 5 knots and 10°, 
from observed winds wherever possible, but supplemented 
by geostrophic winds when necessary. 

The mean temperature at a point, such as point 5 in 
figure 1, is defined as: 
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Ficure 1.—Grid employed in vertical velocity computations, 


The mean advection at the same point is defined by: 


where 


Cus 


+ Ty) + 


Here uv, and v, are, respectively, the west-east and south 
north wind components corresponding to the mean veetor 
wind in box A. C is a proportionality factor which in- 
volves the grid size and the conversion of units. 

Averages over time (12 hours) of advection and tem 
perature change were accomplished by averaging the aé- 
vection fields on two successive maps and subtracting the 
temperature fields. In other words, if the times of sue- 
cessive maps are ¢, and ¢,+12, 


and 


3. SOME EXAMPLES OF THE RESULTS 


Computations of advection and local temperature 
change in the lower stratosphere (200, 150, and 100 mb.) 
were made for two winter periods, December 17-23, 190, 
and January 1-10, 1956. The latter period had beer 
studied by Brockmeyer [1] by means of the trajectory 
technique. 

Figures 2-5 are examples of some of the results of the 
study. The contour patterns at 150 mb. for December 
17, 1953 and for January 4, 1956 are shown in figure 
2 and 8 along with the corresponding patterns of vertical 
velocity. Figures 4 and 5 contain latitudinal cross 8 
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DEC 53 \ 
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Fieure 2.—Contours and vertical velocities at 150 mb. on 
December 17, 1953. 


tions on the same dates. In all four figures, the vertical 
velocities are averages of two successive sets of computa- 
tions; i.e., they are 24-hour averages. 

On December 17 (fig. 2) there was a pronounced ridge 
near 115° W. longitude and a trough near 75° W. Be- 
tween these two meridians was a zone of positive vertical 
velocities above 200 mb. The magnitude of the vertical 
velocities was not large but a smooth pattern is indicated 
and it persisted for the next several periods. The evi- 
dence from the cross section (fig. 4) is that there was 
slight horizontal convergence in this region. Although 
the vertical velocities were still positive (or zero) even 
at 200 mb., the pattern suggests that within the tropo- 
sphere the air was probably subsiding as it moved from 
ridge to trough. Thus, it appears there was a reversal of 
the sign of the vertical velocity between the troposphere 
and the lower stratosphere. 

The case of January 4 (fig. 3) tells a different story. 
Here there was a ridge near 105° W. longitude and a 
trough near 80° W. In this case, however, the region 
between these meridians was one of pronounced sub- 
sidence. The cross section (fig. 5) shows that the sub- 
sidence increased downward. There was no change of 
sign at the tropopause. 

Meteorologists have learned to expect subsidence and 
convergence in the upper troposphere in association with 
flow of air from ridge to trough. Neither of the cases 
presented here, in spite of their differences, conflicts with 
this concept. In both cases, the indicated sign of the con- 
vergence in the lower stratosphere is in agreement with 
that which we would expect in the upper troposphere. 
In one case this convergence resulted in decreasing sub- 
sidence with increasing height; in the other the vertical 
gradient of vertical velocity (0w/dz) appears to have 
been of such a magnitude that there was a reversal of the 
sign of the vertical velocity between the upper troposphere 


Figure 3.—Contours and vertical velocities at 150 mb. on 
January 4, 1956. 


and the lower stratosphere, and the rising motion in the 
lower stratosphere increased with height. 

We are confronted, then, with one case in which the 
stratospheric pattern of vertical velocities was that with 
which we are familiar in the troposphere ; i.e., subsidence 
in air moving with a northerly component; and another 
in which the sign of the vertical velocity was reversed. It 
should be pointed out here that the two examples pre- 
sented here are not extreme cases, although they were 
chosen to illustrate the wide range of vertical velocity 
patterns. It should also be pointed out that all possible 
combinations of vertical velocity and meridional flow were 
encountered, particularly during the period of December 
1953. 


4. THE LONG-WAVE PATTERN 


In an attempt to create some understanding out of the 
seeming disorder apparent on the daily maps, each se- 
quence was separated into a time-mean component and 
deviations from the mean. The means were computed by 
averaging the values at each grid point, and here are 
identified with the long-wave systems which are indeed 
very nearly stationary. The deviations from the means 
are identified with the short-wave systems. In this sec- 
tion we shall discuss only the long-wave (mean) 
components. 

The map in figure 6 illustrates the mean contour, iso- 
therm, and advection fields at 150 mb. for the period 
December 17-23, 1953. (The mean advection is determined 
not by V- 77, but rather as the average of all the fields 
of advection (V: VT). Any correlation between the devia- 
tions of V and V7’ from their means will cause differences 
between these two types of averages. Thus, the field of 
advection in figure 6 (and also in fig. 7) need not agree 
with that which would be given by the contours and iso- 
therms on the maps). Over this period there was a slight 
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decrease in temperature which occurred almost uniformly fore not included in the determination of the mean ver pe 


over the entire area. This decrease in temperature is tical velocity. One then writes: 
taken to be a diabatic radiational heat loss and is there- w= —0.24 V-97 (3) wi 
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Fieure 6.—Mean 150-mb. map for the period December 17-23, 
1953. The isopleths of mean temperature advection are labeled 
in °K. (12 hr.)* and also as mean vertical velocities in em. 
sec.“ according to equation (3). The underlined values are the 
vertical velocities. 


Certainly, if the long waves are indeed stationary and in 
a steady state, there should be no local temperature 
changes which could be associated with them. Thus, the 
isopleths of advection in figure 6 are labeled at one end in 
units of °K. (12 hr.) and at the other in units of em.sec."*. 

The mean map for the period January 1-9, 1956 is 
given in figure 7. In this case, it should be pointed out, 
the averages of the local temperature changes, although 
everywhere small (mostly less than 14° K. per 12 hours), 
were not uniformly negative and did exhibit some sem- 
blance of a pattern. This is due to a combination of 
effects: the fact that an integral number of short waves 
need not have passed a given location during the period; 
non-uniform radiational effects caused by differences in 
the temperature and nature of the earth’s surface and of 
the air itself; also, a slight change of the long-wave posi- 
tion may have occurred. Nevertheless, the local tempera- 
ture changes were sufficiently small that they are neglected 
in figure 7, and equation (3) was employed directly as in 
figure 6. 

The fact that the average local temperature changes 
were so small and uniform offers a considerable degree 
of “a posteriori” justification for identifying the mean 
fields with the atmospheric long waves. If the long waves 
had moved significantly during the periods studied they 
would have carried with them their temperature fields, 
ind this movement would have shown up in the average 
temperature derivatives. Examination of the daily maps 
of the deviations of the contour fields from their means 
(not shown here) also supports the contention that the 
long waves were very nearly stationary during these 
periods. 

Both situations exhibit certain fundamental features 
Which we have come to recognize as characteristic of the 


Figure 7.—Mean 150-mb. map for the period January 1-10, 1956. 
Isopleths are labeled as in figure 6. 


lower stratosphere outside the polar vortex, namely, warm 
troughs and cold ridges, with the isotherms showing 
greater amplitude than the streamlines. This corre- 
sponds, as Kochanski pointed out, to air descending from 
ridge to trough and ascending from trough to ridge. 
Kochanski’s results indicated an average vertical velocity 
of 1.4 em. sec. near the point of inflection at 140 mb., a 
value which is in reasonable agreement with those shown 
in figures 6 and 7. 


5. SHORT-WAVE COMPONENTS 


One may obtain an array of deviations from their 
means of VeV 7,0 7'/dt,w, and h (the height of the 150-mb. 
surface) by subtracting from the daily mean values the 
appropriate means as shown in figures 6 and 7 (and also 
the mean values of 07'/0¢). These deviations are taken as 
representing the effects on each of the parameters of the 
migratory short-wave disturbances. To study the inter- 
relations of the various terms, correlation coefficients were 
computed. The results are summarized in tables 1 and 2. 
The variables correlated are shown by the subscripts to r 
according to the following scheme: 1=VeV7', 2=07'/0t, 
3=wi4=A. 

The largest coefficients are clearly those between ad- 
vection and vertical velocity. The positive value means 
here that cold air sinks and warm air rises; i.e., a direct 
circulation. The principal significance of these coeffi- 
cients is that the local temperature change, which is also 
strongly correlated with the advection, is not closely re- 
lated to the vertical motion. 

On the basis of the correlation coefficients it becomes 
possible to suggest a pattern of the relationships between 
the short-wave troughs and ridges and the corresponding 
components of the vertical velocity. In the vicinity of 
the short-wave troughs there tends to be a maximum of 
cold air advection and local cooling. The advection and 
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TABLE 1.—Correlation coefficients computed for the period December 
17-22, 1953. Subscripts on r refer to deviations from their mean 
values of the variables as follows: 1=V-V7 T, 2=0T/0dt, 3=w, and 
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TABLE 2.—Correlation coefficients computed for the period January 
1-9, 1956. Subscripts on r refer to deviations from their meq, 
values of the variables as follows: 1=VV-T, 2=0T/Olt, 3=w, and 


4h. nis the number of points on which the respective coefficients 4=h. nis the number of points on which the respective coe ficients 
are based. are based 
n ra ria ro ro n T13 123 | Ty 
An 81 |—0.727 |—0.906 |—0.537 | 0.368 |-0.500 9, 
40° 442 |-0.816 |—0.416 |—0.158 | 0.158 | 0.170 —0. 906 |—0. |—0. 
60 | —. 49 | —.509 |) 556 367 99 | 682 | —.908 | —.324 | | —.380| 
66 | —.761 | —.767| —.676| .167| .749| 90 | | —. 850 | —. 285 | | —.276 | ay 
66 | —. 767 | —.446 | —.739 | — 101 72| —.715 | —.850 | —.549 | —.309 
All latitudes: All latitudes: 2 
dee. 41 | —.502 —.920| —.916| .146| .486| Jan. 38 | —.662 | —.957 | —.331 .414) 
Dee. 18 41 | —.040 | —. 733 | —.020 | —.651 | .645| —.424 Jan. 38 | —. 763 | —.943 | —.435 | | 
Dee. 19-- 41 | —.576 | —.852| —.112| .062| .244| —.020 Jan. 3__- 38 | —.839 | —.833 | —.609 | .397 | 673 | 
Dec. 20 41 | —.387 | —.704 | —.422| —.254 .365| Jan. 4... 38 | —. 684 | —.583 | —.199 | —.194 | 
Dee. 21. 41 | 768 | —. 569 | —.823 | -.097 . 655 434 Jan. 5. 38 | —.158 | —.905 | —.377 | —.277 | | 
Dec. 22 | 41 | —.696 | —.848 | —.527| .200| .755 159 = re 
All days and all latitudes 246 —.616 | —.7 —. 502 | —. 036 592 267 38 | — 796 — 678 303 ‘004 080 
Jan, 88 | —. 165 | —.926 | —.550] —.218 | 
All days and all latitudes__..| 342 | —.684 | —.882 | —.394 259 353 
the local temperature change, being of opposite sign, 


tend to cancel, and the vertical velocity associated with 
the short wave is consequently of small magnitude. 
However, since the magnitude of the advection in gen- 
eral exceeds that of the local temperature change, there 
is a tendency for the air to sink as it passes through 
troughs. Whether the vertical velocity which is here 
associated with the migratory disturbances is sufficiently 
persistent, or of great enough magnitude, to be of synop- 
tic significance is doubtful. However, in spite of the 
small magnitude of the correlation coefficients, 7,, the 
agreement between the two independent periods lends 
credence to the pattern presented above. 

In principle, variables which are poorly correlated may 
still be highly dependent upon one another, if the vari- 
ables are 90° out of phase. Additional study of the data, 
however, indicates that this is not the case with regard 
to the relationship between vertical velocity and contour 
height. 

6. SUMMARY 


The vertical velocities in the lower stratosphere do not 
present the consistent pattern which is normally found 
in the upper troposphere. Flow from ridge to trough 
_sometimes has a positive and at other times a negative 
component of vertical velocity. Nevertheless, the pat- 
terns of horizontal convergence in the lower stratosphere 
appear to be continuous with those in the upper 
troposphere. 

The stationary and transient waves have been treated 
separately. The present evidence supports the findings 
of Kochanski [3] that, in association with the long- 
wavelength, quasi-stationary waves, the air subsides in 
flowing from cold ridges to warm troughs. Study of the 
shorter-wavelength, transient waves indicates a consid- 
erable degree of correlation between temperature advec- 
tion and the height anomaly (—0.592 and —0.394 in the 
two independent cases), between local temperature 
change and height anomaly (+0.592 and +0.353), and 
between advection and local temperature change 


(—0.616 and —0.684). There is also a lesser correlation 
between vertical velocity and height anomaly (+0.267 
and +0.293). 


These results suggest a model in which 


there is considerable cold advection and local cooling, 
and a lesser degree of subsidence, in association with the 
short-wave troughs. 

By superimposing the models suggested above, one ar- 
rives at a picture of the flow in the lower stratospher 
which is in reasonable accord with observations. In par- 
ticular, it is found that all possible combinations of posi- 
tive and negative advection, local temperature change, 
vertical velocity, and meridional flow are to be expected, 
depending on the relative positions (and presumably als 
the intensities) of the long- and short-wave disturbances 
Nevertheless the vertical velocity will always be well cor- 
related with the advection; that is to say, the effect of 
local temperature change is small, by comparison, in the 
short waves, and nonexistent (or almost so) in the long 
waves. For the short waves the advection and vertical 
velocity were correlated by —0.765 and —0.882. For 
many purposes, then, the temperature advection, whieh 
can be determined from a single map, can be used as an 
adequate measure of vertical velocity. 
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THE RELATION OF WEATHER FACTORS TO THE YIELD 
OF WINTER WHEAT IN BOX ELDER COUNTY, UTAH 


MERLE J. BROWN 


Weather Bureau Airport Station, Salt Lake City, Utah 


(Manuscript veceived December 29, 1958; revised February 16, 1959] 
ABSTRACT 


A multiple linear regression equation is developed which relates the yield of winter wheat to “effective” pre- 
cipitation during certain months of a crop year. Antecedent precipitation and other weather factors are in- 


vestigated to determine their effect on yields. 


1. INTRODUCTION 


The relation of weather factors to crop yields is be- 
ing studied now more than ever before. The economist 
is interested in learning more about these relationships so 
that future economic conditions can be foretold more 
accurately. Agronomists are concerned with such studies 
so that, among other things, wiser recommendations for 
crop planting can be made for various areas. The farmer, 
the geographer, the crop insurance worker, and others 
also benefit as our knowledge of weather-crop relation- 
ships increases. In this paper an attempt is made to cor- 
relate the yield of one of the important dryland crops in 
Utah, winter wheat, with certain weather factors. 

It is a well-known fact that wheat is the most important 
grain crop in the world. Wheat is adapted to a wide 
range of climate and soils, but it is grown principally 
inmoderately dry temperate climates having am average 
amual precipitation of less than 30 inches [4]. The 
highest average yields are usually obtained in regions 
with an annual precipitation of 25 to 30 inches [3]; but 
in Utah the principal dryland farming areas receive less 
than 15 inches, on the average [9]. 

There are several weather factors affecting the growth 
and yield of wheat. Zink [11], in his study of the weather 
elements important to wheat growth and yield near Nephi, 
Utah, found the highest correlations were with evapora- 
tion, length of drought periods, and rainfall. He also 
found that correlations with mean temperatures were low, 
but were higher with maximum and minimum readings. 
Zink concluded that the best statistical results are obtained 
by using plant growth stages. However in the present 
study which is confined to Box Elder County, Utah, little 
ifany phenological data for winter wheat were available, 
and it was necessary to use weather data for specific calen- 
(ar periods. Also, there were no evaporation records 
aailable for use in this study. 

Some of the factors that other investigators have be- 
lieved to be important are: the amount, distribution, 


reliability, and effectiveness of rainfall; length of grow- 
ing season; the amount of sunlight and soil moisture; 
snowfall; snow cover; and the frequency and duration 
of very low or very high temperatures [1, 11]. 

The principal climatic factor limiting the yield of 
winter wheat over western United States is insufficient 
precipitation. In the growth of this crop there appears 
to be at least two critical periods in which adequate soil 
moisture is essential to the production of maximum 
yields. One of these is the period during heading, bloom- 
ing, and filling. Another critical period is the time of 
planting [8]. 

Most of the winter wheat in Utah is of the hard red 
class and is grown principally on “summer fallow” 
ground. This crop is usually planted in September and 
harvested in July. Diseases such as rust, scab, mildew, 
and leaf spots; leaching of plant food from the soil; ex- 
cessive growth and lodging; and difficulties in planting, 
harvesting, and caring for the crop are not as much of a 
problem in Utah as in more humid climates [6]. 

Although this particular study is confined to Box Elder 
County, which ranks first in the production of small 
grains in Utah, it is believed that the results would be 
applicable to other dryland winter wheat areas of western 
United States. Box Elder County, in northwestern 
Utah, comprises an area of 5,600 square miles, with about 
90,000 acres planted to dryland winter wheat [7]. Most 
of the cultivated crops are produced in the southeastern 
portion of the county. The northern and western sec- 
tions are used principally for grazing. 


2. PROCEDURE AND RESULTS OF STUDY 


Winter wheat yield data were obtained from the offive 
of the Agricultural Marketing Service, U.S. Department 
of Agriculture, Salt Lake City, Utah; temperature and 
precipitation data were taken from Climatological Data, 
Utah (monthly and annual issues) published by U.S. 
Weather Bureau. 

The first step was to analyse the yield data for time 


ry 
In 4 
ud 
153 
024 
‘ 
25; 
301 
Oe 
42 
he 
re 
1 
) j 
| 
of 
he 
ng 
ral 
‘or 
ch 
all 
nd 
he 
yas 
nd 
ato- 
Jni- 
mb 
tor, 
Air 
the 
yor! 
ice, 
ver 
vol. 
4 
| 


98 MONTHLY WEATHER REVIEW 


WINTER WHEAT YIELD IN BUSHELS PER ACRE 


1925 ‘x0 ‘40 ‘a5 ‘50 1955 
YEAR 


Ficure 1.—Winter wheat yields, in bushels per acre, for Box Elder 
County, Utah, plotted against years for the period 1926-55. 
The regression line is dashed. 


trends during the period under consideration (1926— 
1955) to determine the effect of improved cultural prac- 
tices, etc., on average yields. According to a study made 
in 1951 by the U.S. Department of Agriculture [10], the 
average yield of a// wheat in the United States showed 
a pronounced increase during the period 1939-1948, com- 
pared to the eight decades prior to 1939. In the Western 
States, an early peak was reached in 1915, with the yield 
dropping markedly in 1917 and remaining relatively low 
through 1940. 

In Box Elder County there is no evidence of a marked 
time trend in the winter wheat yield data from 1926 
through 1955. A linear regression line fitted to yields 
shows a slight decrease of about 1 bushel per acre during 
this period, but a “t” test indicates this change is not 
significant at the 5 percent level. The regression line 
and a plot of annual yields versus years are shown in 
figure 1. 

Precipitation records for ten weather stations in and 
near Box Elder County were tested for consistency for 
the period 1926-1955, by the double-mass analysis tech- 
nique [2]. Of the ten stations, seven were found to have 
compatible records that could be used without any adjust- 
ment of the data. The correlations between seasonal pre- 
cipitation (July of one year through June of the follow- 
ing year) at five of the stations (Grace, Idaho; Brigham 
City, Corinne, Richmond, and Riverdale, Utah) and the 
annual wheat yields in bushels per acre for Box Elder 
County were computed. From examination of these cor- 
relations and some analyses by multiple regression, Grace 
and Riverdale were dropped and precipitation at Brigham 
City was assigned a weight of 0.50, Corinne 0.30, and 
Richmond 0.20. 

The effective monthly precipitation for the period under 
consideration was then computed. Effective precipitation 
is the sum of the precipitation values times the respective 
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TaR_e 1.—Predicted and actual yields of wheat in bushels per acy 
for Box Elder County, Utah, for 1956 and 1957 


Year Predicted | Actual | /redieted 
yield yield ~ actual 
yield 


station weights. This procedure is identical with the one 
outlined by Linsley, et al. [2], except effective precipita. 
tion in this study is correlated with wheat yield rather 
than water-year runoff. 

Since it was known that for the production of maximuy 
winter wheat yields, adequate soil moisture is necessary 
at the time of planting and during heading, blooming 
and filling, precipitation data during the months of Sep. 
tember, October, May, and June were used in this study, 
A three-variable correlation was computed, using as pa- 
rameters the following: Y=annual yield of dryland win. 
ter wheat in bushels per acre; Y,=effective precipitation 
in inches for September plus October; and X.=effectiye 
precipitation for May plus June. The relationship of 
these variables is shown by the scatter diagram of figure?. 
The sloping lines on this chart are least squares lines for 
annual winter wheat yield. The regression equation was 
determined to be : 

Y=8.8+1.67 X¥,+1.85 

where Y=estimated annual yield of winter wheat. The 
coefficient of multiple correlation # for this set of data 
is 0.735. # measures the degree of linear association be- 
tween the estimated and observed values of the dependent 
variable [5], in this case between the estimated and ob- 
served average wheat yields. *, that part of the total 
sum of squares which has been explained by linear re 
gression, is 0.540. In other words, about 54 percent of 
the variation in wheat yield can be accounted for by var- 
ations in. XY, and XY». 

The average forecast error for the sample was 2) 
bushels per acre. Yields were forecast too low 17 times 
with an average error of 2.2 bu. per acre, and were fore 
cast too high 13 times with an average error of 2.5 bu 
per acre. The standard error of estimate was computed 
as 2.73. 

Average wheat yields for Box Elder County, Utal, 
are not available prior to 1926. 30 years of record wis 
considered a minimum amount of data for a study sue 
as this; hence, only 2 years (1956 and 1957) of test data 
were available. A comparison of the predicted and actu 
wheat yields for those years is shown in table 1. 


3. OTHER FACTORS INVESTIGATED 


Winter kill of wheat is occasionally a problem in Utal 
[9] and, obviously, this reduces the average yield in some 
years. The amount of damage to winter wheat as a rest! 
of low temperatures depends, in part, on the actual values 
of the minimum temperatures observed and whether there 
is a snow cover on the ground when the low temperature 
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occur. However, the introduction of these variables into 
the reg ession equation did not improve the results signifi- 
cantly. (The addition of the lowest winter temperatures 
at Corinne, Utah, accounted for only an additional 2 
percent in the variation of winter wheat yields.) 

As was mentioned above, most of the winter wheat in 
Box Elder County is planted on “summer fallow” ground. 
One might expect that the inclusion of an antecedent 
precipiation variable in the regression equation would 


improve the predictability of winter wheat. However, 


the addition of the following variables did not change the 
results significantly : (1) The sum of the effective precipi- 
tation amounts for April, May and June of the previous 
cropyear. (2) Effective precipitation for the entire previ- 
ous crop year. 

This is not to say that summer fallowing is not effective. 
It may be that enough moisture is added to the soil, even 
in the driest years, so that the amount of precipitation 
during the fallow year is seldom the principal limiting 
factor in the production of maximum yields, or possibly 
the measure of “fallow year” precipitation used does not 
adequately express the important factor of soil moisture 
at planting time. 

The relationship of other weather factors to the yield 
of winter wheat was investigated. Some of these were: 
(1) The lowest monthly winter mean temperature at 
Corinne, Utah. (2) The mean temperature for the months 
of April, May, June, and July at Corinne. (3) The effec- 
tive precipitation for the entire crop year, August through 
July. (4) Maximum temperature for the months of 
May and June at Corinne. 


4. CONCLUSIONS 


Dryland winter wheat yield in Box Elder County, Utah, 
is related principally to the amount and seasonal distri- 
hution of precipitation. 

The results of this study indicate the importance of 
adequate moisture in September, October, May, and June. 
If sufficient soil moisture is available during those months, 
there is an excellent chance of an above average yield in 
the crop year. 

The addition of other variables to the regression equa- 
tion, such as effective precipitation during the previous 
crop year, winter minimum temperatures, and effective 
precipitation for the current crop year, did not improve 
the results significantly. 
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EFFECTIVE PRECIPITATION FOR MAY 


L 
© 2 SO 75 100 25 SO 75 200 25 SO 75 300 25 SS 75 400 25 S 75 500 
EFFECTIVE PRECIPITATION FOR SEPTEMBER + OCTOBER (INCHES) 


Ficure 2.—Joint relationship between September plus October 
precipitation and May plus June precipitation and yield of 
winter wheat, in bushels per acre, Box Elder County, Utah. 
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Weather Note 


SNOWFALL AT SAN ANTONIO ON MARCH 16, 1959 


HAROLD S. McCRABB 


Weather Bureau Airport Station, San Antonio, Tex. 


[Manuscript received April 13, 1959] 


Snowfall at San Antonio in any month of the year is an unusual 
occurrence—an event for the local populace. A flake or two of 
snow falling in San Antonio causes the telephone lines to be 
jammed at the Weather Bureau, is headlined in the newspapers, 
and becomes the principal topic of conversation. The total num- 
ber of days of measurable snowfall at San Antonio from 1885 to 
1958, inclusive, is 27, with an average fall of about one-half inch 
a year. The average for March is a trace. 

On the morning of March 16, 1959 hundreds of people called the 
San Antonio Weather Bureau office to report the falling of snow. 
A quick look at the meteorological conditions normally would have 
suggested discounting any reports of snow, as the current tem- 
perature was 50° F., and the freezing level was about 6,000 feet 
above the surface. Yet from 0755 to 0830 cst on March 16 flakes 
of snow were observed and recorded as a trace at the Weather 
Bureau and corroborated by a number of experienced observers. 
The ceiling at the time of the snowfall was 9,000 feet, the lowest 
reported at any time that day. San Antonio raobs indicated the 
freezing level at 8,000 ft. mst at 0548 cst and at 6,200 ft. Mst at 
1737 cst. A perusal of all available data indicated that the freez- 
ing level was at least 6,000 feet above the surface at the time of 
the snowfall and the surface temperature never fell below 50° F. 

Conditions aloft were favorable for the occurrence of precipi- 
tation. At 700 mb. a trough passed over the station. At 500 mb. 
there was evidence of “digging” to the southwest of San Antonio 
with a major trough shearing over North Texas, and the forecast 
cut-off Low in West Texas was evident on the 1200 emt 500-mb. 
chart for March 17. The temperature at 500 mb. dropped from 


—16° C. to —19° C. in the 24 hours ending at 0000 emt, March 17. 
The 250-mb. chart for 1200 emt, March 16 showed two SW-NE jet 
streams, one over North Texas and one over South Texas. The 
chart 12 hours later showed one jet stream through South Texas, 


and later charts showed a continued southward shift of the je 
stream. A normal sequence of events followed—a wave was ip- 
duced and developed on the cold front in the Gulf of Mexico. 

From thermodynamical considerations formation of snow aloft 
in this case was logical, but for it to reach the ground under the 
conditions on the morning of March 16 does not appear to be s 
reasonable. It is a common occurrence for aircraft to encounter 
snow in clouds above the freezing level, and for no snow to reach 
a surface with temperature well above freezing and far below the 
freezing level. When snow falling through the atmosphere melts, 
it absorbs from the air the same amount of heat as was released 
when it was being crystallized. Certainly, with very little snow 
falling through a thick layer of air with temperatures well above 
freezing, there should be no difficulty in the snow melting before 
it reaches the ground. In fact, it would seem impossible for snow 
under such conditions to reach the ground. 

The most plausible explanation of the March 16 snowfall at San 
Antonio appears to be that the snow reached the ground in the 
cool air of vertical downdrafts from high-level thunderstorms or 
cumulonimbus clouds. There was some evidence of this from 2 
slight increase in the surface winds and an increase in the sur- 
face pressure at the time of snowfall. Pilot reports indicated 
that a high-level instability line moved concurrently with the 
upper trough. Radar and pilot reports indicated tops of clonds 
to 31,000 feet. There were a number of pilot reports of snow, 
sleet, freezing rain, and hail at 4,000 to 8,000 feet ms east of San 
Antonio as the trough and high-level line of thunderheads 
thunderstorms continued their eastward movement. Beeville re 
ported sleet and rain at 0900 cst; Victoria reported sleet at 001 
cst; and Houston reported small hail at 1010 cst. Before the 
snowfall at San Antonio and before sunrise, snow was reported in 
West Texas and eastern New Mexico. 
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ABSTRACT 


Most meteorological soil moisture budgets do not account for soil moisture stress changes in the drying 
cycle or for changes in ground cover or expanding root system. A simple modulated technique is described 
which considers these factors. Soil moisture stress determined by this technique has a significantly higher 
correlation with wheat yield than does moisture stress determined by a common method. 


1. INTRODUCTION 


An important phase of research in agricultural meteor- 
ology and soil physics is concerned with soil moisture and 
the process of evaporation. The soil, plant, and at- 
mosphere each play important roles in this process but 
are only parts of a continuous system of transport of 
water from the earth to the atmosphere. While much 
emphasis has been placed on the individual components 
of this continuum, there has been less attention devoted 
to the system as a whole than it would appear to merit. 

It is recognized that moisture use by crops cannot pass 
a maximum which is determined by the amount of energy 
available, from the sun and the horizontal transfer of 
heat, to convert liquid water to water vapor. As long 
as water is freely available for evaporation at the earth- 
atmosphere interface, the rate of evaporation is mostly 
dependent upon meteorological parameters. The amount 
of evapotranspiration required to satisfy the atmospheric 
demand has been called potential evapotranspiration 
(PE) [20]. However, as the soil dries, the available 
moisture decreases, hydraulic tension increases, and trans- 
port of water to the interface is insufficient to meet the 
atmospheric demand and actual evapotranspiration (AE) 
falls short of PE [6,7]. 

The commonest soil moisture budget consists of sub- 
tracting daily PE from daily rainfall. This amount is 
then subtracted (or added when rain exceeds PE) from 
the moisture present in the soil to give the new soil mois- 
lure storage, until available stores have been exhausted. 
For budgeting irrigation water, this process is continued 
ona daily basis until a certain permissible soil moisture 
deficit occurs. At this time an amount of water equal to 
the deficit: is applied to return soil moisture storage to 
lield capacity. This isa simplified description of moisture 
sé by crops under irrigation. “Dryland” soils present 


many problems with meteorological moisture budgets be- 
cause they are rarely at field capacity and plant roots 
explore deep moisture stores. 

Most methods of determining PE are semi-empirical 
[2, 6, 10, 14, 15, 20,21]. Evaporation pans or atmometers 
have also been used [3, 5, 8, 12,24]. As stated previously, 
when the soil dries out AE becomes less than PE, 
Thornthwaite suggests that AE is in ratio to the soil mois- 
ture in storage (SM). That is, when SM is \% the total 
storage possible, AE is 4% PE. He proposes that evapo- 
transpiration continues to near oven dryness. Thorn- 
thwaite’s data supporting his soil moisture depletion 
curve were obtained from vapor pressure and tempera- 
ture profile measurements taken at O'Neil, Nebr. [21]. 
Curves C and D, figure 1, represent Thornthwaite’s con- 
cept; curve D continues to the wilting point and C to oven 
dryness. 

Blaney and Criddle [2] correlated actual measure- 
ments of consumptive use (evapotranspiration) with 
monthly mean temperature and daylength in an attempt 
to obtain a soil-plant coefficient that, when incorporated 
into their empirical relationship, would give it a wide 
application. The plant coefficients varied from 0.85 for 
alfalfa to 1.2 for rice, and applied to semiarid and arid 
conditions. : 

Pierce [16] used Thornthwaite’s method of calculating 
PE and compared these data to total water use by 2d 
year meadow as measured by weighing lysimeters. 
Thornthwaite’s PE was adjusted upward to correct for 
crop density, age, and rooting depth. Pierce’s “dryness” 
correction curve is represented in curve B figure 1, while 
his seasonal correction for 2d year meadow (state of 
growth and time of year) is represented in figure 2. This 
curve is thought to apply to many soils, provided drain- 
age is not restricted. 
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Figure 1.—Various proposals for adjustment of PE as soil moisture 
decreases. A, Veihmeyer; B, Pierce; C and D, Thornthwaite. 


Van Bavel’s [22] method of determining PE is a sim- 
plification of Penman’s equation. Nomograms have been 
drawn to simplify the estimation of daily PE. No allow- 
ance is made for soil dryness and crop density (ground 
coverage) and only moisture within the root zone is con- 
sidered. (This amount would increase as the crop de- 
velops with increasing depth of root zone.) Curve A 
figure 1 represents Van Bavel’s relationship between soil 
dryness and PE. This curve is similar to that proposed 
by Veihmeyer [24]. 

In a recent review of Russian literature, Lemon [11] 
presents several soil moisture evaporation curves. They 
are divided into three portions; (2) evaporation proceeds 
in accordance with the atmospheric demand, (b) evap- 
oration rate declines rapidly as moisture films to the 
surface become discontinuous and transfer of moisture 
to the interface decreases, (c) extremely slow moisture 
movement is dominated by adsorptive forces at liquid- 
solid interfaces within the soil. Curves in figure 3 sug- 
gest this concept. Except for the initial plateau, this is 
characteristic of many “tension-moisture content” curves. 

Marlatt [13] investigated the change in AE as the soil 
in lysimeters and field plots dried out. By regular soil 
sampling at 3-inch intervals to 48 inches under a corn 
crop throughout the season, curves similar to those of 
Lemon [11] were obtained. He found that AE pro- 
ceeded at the potential rate up to a point depending 
chiefly on rooting depth, then AE fell off sharply. These 
data are represented by the curves in figure 3. The de- 
flection points correspond to the various crop rooting 
depths at different periods during the season. Philip 


[17] obtained the same type of curves by analysis of the 
mathematical and physical aspects of soil moisture evap- 
oration. The horizontal portion of the curves may 
roughly correspond to the zone of “complete depletion,” 
described by Hagan [4]. 


This zone may be defined as 


100 + 
First Cutting 
Second Cutting 
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Ficure 2.—Adjustment of PE for stage of growth of meadoy. 
(Pierce) 


the volume or, in the case of closely spaced plants, the 
depth from which soil moisture is withdrawn at or near 
the potential rate until most of the available water has 
been removed. The soil in this zone would be thoroughly 
permeated with fine roots. The second and third portions 
of the curves may correspond to the “transitional zone” 
of Hagan [4]. This area contains a varying concentra. 
tion of roots, decreasing with distance from the plant 
and extending to a considerable depth. 

Marlatt [13] compared the measured and computed 
soil moisture content of lysimeters and _ field plots 
Measurement was accomplished by sampling, and com- 
putation was achieved by adjusting Thornthwaite’s PE 
to soil dryness and rooting depth (according to curves 
in fig. 3). His correlation was 0.997. Using Van Bavel’ 
[22] nomograms he found departures of the computed 
and measured SM, particularly when the soil was dry. 


2. MOISTURE BUDGETS 


To be useful, a scheme for the budgeting or control of 
soil moisture by the meteorological method should have 
several characteristics. First, the method of characteriz- 
ing the drying ability of the atmosphere should (a) be 
convenient and simple to use, (b) integrate into one 
measurement all the various meteorological factors 
affecting the evaporating ability of the atmosphere, (¢) 
(instruments) be free of structural components that dis 
tort evaporation measurements. 

Second, the manipulation of evaporation data shouli 
(a) be realistic in its description of natural processes, 
(b) have accuracy compatible with the use to which it 
is put, (c) be simple and practical. 

Many of the methods presently used to measure evap 
oration do not fill the above requirements as they are nol 
simple or convenient. Those that are, do not accurately 
integrate meteorological factors affecting evaporation inte 
one measurement. Many instruments are cumbersome and 
are not free of structural disadvantages. A case in poill 
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Fievre 3.—Adjustment of PE for soil dryness and rooting depth of 
crop (Lemon, Marlatt). Curves A to D correspond to increases 
in rooting depth of crop. 


is that of evaporation pans. Many authors have claimed 
that atmometers fill most of the above requirements [5, 
12, 24], and two recent papers compare an evaporation 
pan and several atmometers [8, 18]. The black Bellani 
plate atmometer was described as an accurate and simple 
instrument that adequately characterizes the evaporation 
demand of the atmosphere for water. 

Similarly, few of our present moisture budget methods 
(as opposed to evaporation measuring techniques) fulfill 
all of the above requirements. The accuracy of many 
present techniques is doubtful; most of them are not 
realistic in that they do not account for the drying of the 
soil as moisture loss continues. 

It is clear that any accurate description of the moisture 
status of a soil will not be simple. The plant is the only 
true indicator of this factor and at the present ‘time it is 
not possible to measure plant moisture stress, per se. It 
isnecessary then in any budget scheme to make some sim- 
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plifying assumptions, which will be based on crop growth 
stage, rooting characteristics, and other soil and plant 
factors. In making these assumptions it is desirable, at 
the same time, to fulfill the requirements set forth con- 
cerning evaporation data manipulation. 

A soil moisture budget has been developed which, in the 
opinion of the authors, fulfills the requirements of the 
“model” budget described above. Soil dryness is com- 
pensated for in a manner similar to that described by 
Lemon [11], Marlatt [13] and Philip [17]. Crop rooting 
depth and soil storage may be adjusted periodically dur- 
ing the growing cycle, and rainfall is evapotranspired at 
the potential rate on the day of the rain and on subsequent 
days depending on the amount of rain. Percolation, run- 
off, and soil dryness by zones are also estimated. 

The assumptions made in the technique are: (a) all 
moisture from the zone ot “complete depletion,” described 
earlier, is evapotranspired at the potential rate, (b) PE 
is modulated so that moisture in the “transitional zone” 
is withdrawn at a decreasing rate depending on the per- 
centage of available moisture remaining, (c) available 
moisture is withdrawn from the topmost layer of soil, 
before extraction occurs from lower layers, (d) rainfall 
occurring after a dry spell is evapotranspired at the po- 
tential rate until it has been depleted, unless rainfall per- 
colates into the “transitional zone”; then PE is modulated 
and withdrawal is at a slower rate, (e) percolation is 
regarded as complete; if the soil is saturated, rainfall is 
regarded as runoff. Soil moisture between saturation and 
field capacity is assumed to percolate. 

The weakest assumption is (c) since a plant withdraws 
water from any zone occupied by roots so long as soil 
moisture in that zone is available. However, the essen- 
tials of the problem are retained since the plant withdraws 
water from soil zones of varying moisture contents at 
varying rates, depending on the moisture content of each 


TaBLe 1.— Daily PE and soil moisture 


Soil moisture storage—Percent of PE 
100 50 20 10 | 5 
Month—April 
Soil zone caparity (inches) 
0. 25 0. 25 0. 25 0. 25 0. 25 
Date PE Rain AE SM AE SM AE SM AE SM AE SM Deficit | Surplus 

cee 2.08 4 25 25 25 25 10 01 24 0.00 25 0.00 0.00 
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I'icure 4.—Soil moisture loss calculated by modulated budget (AE) 
and by simple budget (PE). (A indicates spot Coleman moisture 
block readings. 


zone. Recent work by Vasquez and Taylor [24] substan- 
tiates this view. There may be some criticism of assump- 
tion (e) since local conditions and individual soils may 
seriously retard percolation. 

The individual operations of this technique are set forth 
in table 1. For purposes of demonstration, the soil zone 
of complete depletion was assumed to contain 0.25 inch 
of available water and the transitional zones, 1.0 inch of 
water, making a total soil storage of 1.25 inches of avail- 
able moisture. The uppermost 0.25 inch of moisture was 
evapotranspired at 100 percent of PE; the remaining 1.0 
inch was withdrawn in a stepwise fashion with 0.25 inch 
per step, in such manner that the steps were the best fit 
to the type of curves in figure 3. In other words, the first 
0.25-inch portion of the 1.0 inch in the transitional zone 
was withdrawn at 50 percent, the second at 20 percent, 
the third at 10 percent, and the fourth at 5 percent of PE. 

Figure 4 compares soil moisture loss by a grass-legume 
sward calculated by two techniques: (a) the modulated 
budget described above, and (b) the common type dis- 
cussed briefly at the beginning of this paper. PE was 
estimated with a black Bellani plate atmometer [8]. The 
differences are clearly discernible particularly as the soil 
begins to dry out. Spot moisture block readings verified 
the accuracy of the modulated technique. 


3. COMPARISON OF TWO BUDGET TECHNIQUES 


The two moisture budget techniques partially described 
above were programmed for an IBM 650 and soil moisture 
calculations were made for Lethbridge, Alberta, for 36 
years (1921-56 inclusive). In both techniques, PE was 
estimated by Thornthwaite’s method [19]. 

Common Budget (A): The method employed in the 
simple budget technique has been outlined by Robertson 
and Holmes [18] and others [20, 22]. 

Modulated Budget (B): The basic characteristics of 
this technique are outlined in section 2 of this paper. The 
apportionment of the soil zones and the amount of mois- 
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Figure 5.—IBM 650 program for adjustment of PE as soil dries 
and plant roots expand. (Curve A for Aug. 1 to May 31; curve 
B June 1 to June 30: curve C for July 1 to July 31.) 


ture in each zone is shown in figure 5. The soil drying 
curve at Lethbridge for stubble or trash covered plowed 
ground is approximated in the dashed curve A. The step- 
wise best fit curve used in programming the IBM 650 is 
drawn over the ideal curve. The length of each step 
indicates the amount of moisture leaving the soil at the 
rate of PE shown on the ordinate. This curve was used 
to calculate soil moisture loss from August 1 (approxi- 
mate harvest date) to May 31 (seeding date, approxi- 
mately May 1). The June soil moisture loss curve is rep- 
resented in curve B, and conditions during July are shown 
graphically in curve C. These curves were established 
by laboratory and field observations. Yield data froma 
36-year old “continuous spring wheat” rotation at Leth- 
bridge were available for correlation comparisons. 

Using the common and modulated budgets, PE, AE, 
and soil moisture deficits (monthly average below field 
capacity) during the growth months (May, June, and 
July) were calculated and compared, through multiple 
and simple linear correlation, with yields of wheat. The 
results are shown in table 2. 

The data indicate that wheat yields at Lethbridge wer 
not. significantly correlated with the evaporating ability 
of the atmosphere (PE). Soil moisture deficit was more 
significantly correlated with wheat yields than was actual 
moisture use (AE). Plant response has been shown to 
be more closely related to the energy required to take up 
water (e.g. deficit), than to moisture use per se [4, 13, 17], 
The correlation between yield and average monthly mois 
ture deficit below field capacity improves as the seas0l 
progresses. This indicates that as the crop ages and the 
roots are more fully ramified throughout the soil, such 
factors as rooting habit, soil moisture stress, etc., may 
cause a more pronounced effect on the yield than atmos 
pheric conditions. Allowance is made for these factor 
in the modified budget (B). Army and Ostle [1] noted 
an inverse relationship between the evaporating ability 
of the atmosphere (e.g., PE) and evapotranspiration from 


neo 


} 
4 T 
3 
2 
ON, ak 
’ may JUNE JULY AUGUST SEPTEMBER “he 
-t 
2 Al 
D 
-3 
-4@}- 1 
PE fac 
fe 
a 
| T 
(; 
in 
th 
tu 
of 
in 
et 
4 m 
as 
me 
bu 
ror 
ib] 
crt 
in 
of 
Be 
Cu 
it 
| 
3. 
‘ 4, 
‘ 


Tr & & 


March 1959 MONTHLY WEATHER REVIEW 105 


Tarte 2.—Simple and multiple correlations between soil moisture 
actors AE, PE, and deficit (calculated by common (A) and modu- 
lated (13) methods) and yield of wheat at Lethbridge, Alberta 


Correlation coefficient and method 
Soil moisture May ry: June ry2 July rys Ty 123 
A B A B A B A B 
ppt..----------- —0.27 | —0.27 | —0.04 | —0.04 | —0.14 | —0.14 0. 28 0. 28 
—.37 | *—. 38 32 |**—. 46 | *—. 44 |**—. 63 46 | **.68 


1 PE is calculated the same way in both A and B methods. 

: AE is not calculated in method A. 

ryi,ga ete., indicates correlation between yield and May, June, and July soil moisture 
at the 1 and 5 percent levels respectively. 
wheat (e.g., AE). An explanation of this anomaly was 
found in the complex interrelationship of plant growth, 
available soil moisture, and climate under semi-arid con- 
ditions. Results reported here indicate a similar trend. 
The correlation between August PE and AE is r= —0.57 
(significant at the 5 percent level). 

Even though the results with the modulated budget 
indicate considerable versatility and improvement over 
the simpler meteorological methods of estimating soil mois- 
ture status, one may wonder why correlations with yield 
of wheat are not higher. Hopkins [9] pointed out that 
inhibitory factors such as wind, disease, weeds, insects, 
ete, are important in many seasons, and that statistical 
methods may lead to some underestimation of the actual 
association of yield of wheat and meteorological factors. 


4. SUMMARY 


The results presented show that a realistic, yet simple 
meteorological soil moisture budget is possible. Such a 
budget readily accounts for soil moisture stress and plant 
rooting characteristics, etc., and yet has accuracy_compat- 
ible with other methods of determining moisture use by 
crops. When soil moisture drying curves (such as those 
in fig. 5), and accurate estimates of the drying ability 
of the atmosphere such as those obtained with a black 
Bellani plate atmometer are available, a daily record of 
soil moisture is possible. By changing the soil drying 
curves and coefliients, the scheme can be programmed to 
fit a wide variety of soils, crops, and crop sequences. 
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CORRESPONDENCE 


Discussion of “A Modulated Soil Moisture Budget’’ 


M. A. KOHLER 


U.S. Weather Bureau, Washington, D.C. 


January 26. 1959 


Holmes and Robertson [1] have demonstrated that mean 
monthly moisture deficiency is better correlated with 
wheat yield when it is derived by application of a budget 
which considers moisture stress and variable root zone 
than when it is derived from one which neglects these 
factors. They state that “one might wonder why the 
correlations with yield of wheat are not higher” and cite 
Hopkins’ view that “inhibitory factors such as wind, dis- 
ease, weeds, insects, etc., are important in many seasons.” 
The writer feels, however, that one might question the 
reality of the high correlation found (0.68) in view of 
the assumed linear effect of mean monthly deficiency. Cer- 
tainly the yield would approach zero if the deficiency for 
a brief period of days were sufficiently severe even though 
the monthly and season values were not critical. 

What the authors choose to call a modulated budget is 
similar in important respects to that employed by the 
Weather Bureau in connection with its river forecasting 
program [2,3]. The differences in the two methods are 
largely necessitated by virtue of the fact that we maintain 
the budget for a basin (the order of 100 to 500 sq. mi.), 
while the authors are concerned with a particular plot or 
field. Thus while we use a two-layer system, we have not 
found it necessary to vary the two moisture capacities for 
changes in root zone. 

For the basin, we have and use total runoff (surface and 
percolation) as an item in the budget, whereas the authors 
must assume that runoff occurs only when the soil profile 
is saturated. It should be pointed out that observations 
of rainfall and runoff show that the concept of a satura- 
tion threshold is not sound and results in large errors of 
estimated daily runoff. Using this concept, one is forced 
to select unrealistic values of soil moisture capacity if the 
accounting is to provide the correct evapotranspiration 
over extended periods. 

The authors contend that evaporation pans are “cum- 
bersome and are not free of structural disadvantages” and 
have purported to show in an earlier publication [4] that 
the Bellani plate atmometer provides better results than 
can be obtained from a pan. Their comparative studies 
have been restricted to a 4-foot sunken pan of the type 
used in Canada and the conclusions do not necessarily 
apply to the Class A pan used in the United States which 
is free of heat exchange with the soil. All known evap- 
oration instruments possess decided deficiencies and the 
Bellani plate atmometer is no exception. It can be used 


only during those months when the minimum temperatures 
are above freezing and it requires repeated calibration to 
avoid time trends in the data. 

We have found [5] that free-water evaporation can be 
estimated from that of a Class A pan or from certain 
meteorological factors with about equal reliability, Ip 
either case, our moisture budget computations are based on 
the assumption that free-water evaporation and potential 
evapotranspiration are equivalent. The computation of 
potential evapotranspiration from solar radiation, air and 
dewpoint temperatures, and wind, as well as the moisture 
budget computations, can be performed on an IBM 65) 


computer [6]. 
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Reply 


R. M. HOLMES AND G. W. ROBERTSON 
February 27, 1959 


It has been our experience that while periods of a few 
days of high moisture stress certainly influence the yield 
of spring wheat, the effect is not marked unless the def 
ciency occurs at the time of germination and emergence. 
According to our budget, the Lethbridge plots were never 
depleted of available moisture during the entire 36-yeat 
period, although the supply was low oftentimes. Any 
drought periods that affected yield were of sufficient length 
to be accurately accounted for on a monthly basis. In 
other words, rainless periods of a week or less were over 
come because of moisture held in the low root zone 
However, if the drought persisted for longer periods, evel 
this store of moisture became depleted and the effect 
the yield was marked. 


(Continued on p. 118) 
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A NOMOGRAM TO DETERMINE 
MONTHLY POTENTIAL EVAPOTRANSPIRATION 


T. E. A. VAN HYLCKAMA 


Laboratory of Climatology, Drexel Institute of Technology, Centerton, N.J.* 


{Manuscript received November 20, 1958; revised January 16, 1959] 


ABSTRACT 


A nomogram for the solution of Thornthwaite’s empirical relationships among mean monthly temperature, lati- 
tude, and potential evapotranspiration is constructed. In the form of alignment charts, the nomogram is simple to 
use and has the added advantages of the temperature scale being labeled in both degrees Celsius and degrees Fahren- 
heit and of the potential evapotranspiration scale being labeled in both centimeters and inches. 


1. INTRODUCTION 


Potential evapotranspiration is defined as the amount 
of water that would be evaporated from the ground or be 
transpired by the vegetation if there were no shortage of 
water at any time during the period of observation. The 
concept has proven to be a valuable tool in the classifica- 
tion of climates as well as in the computation of water 
balances [1], [2], [3]. 

Thornthwaite [4] found an empirical relationship 
among the mean monthly temperatures, the latitude of 
the location, and the potential evapotranspiration. The 
determination of the latter involves the use of two tables 
and a nomogram and is rather time consuming, especially 
if series of years for different locations have to be 
computed. 

The nomogram presented here was the result of the 
desire to have a fast, accurate, and convenient method to 
determine the potential evapotranspiration according to 
Thornthwaite’s method. The manuscript was ready for 
publication when the author's attention was drawn to a 
recent article by Palmer and Havens [5], dealing with 
the same subject. The alignment charts described in the 
present paper are easier to use and either alone or in com- 
bination with the graphical technique of Palmer and 
Havens will expedite the computations of potential 
evapotranspiration. 


2. THE THORNTHWAITE FORMULA 


The computation of potential evapotranspiration ac- 
cording to Thornthwaite’s [4] formula consists of three 
steps. First, the heat index (/) has to be found for the 
year and consists of the summation of 12 monthly indices 
(i), where and ¢=the mean temperature of 
themonth in °C. Second, the unadjusted potential evapo- 
transpiration (4’) must be computed according to the 
formula 

E=1.6 (1) 


*Present affiliation: U.S. Geological Survey, Washington 25, D.C. 


in which 
a=0,000000675 7* —0.0000771 0.01792 J+ 0.49239. 
. (2) 


Third, must be multiplied by a factor for daylight (/), 
adjusting / to the mean possible duration of sunlight in 
each month of the year. This factor depends on the lati- 
tude of the location. 

The heat index (J) can be computed from Thorn- 
thwaite’s [4] table 1V, but it is simpler to construct a 
conversion scale where both ° F. and ° C. are inserted. 
Such a scale is presented at the top of figure 1. By 
adding the #’s for each of the twelve months of the year. 
J is determined. This completes the first step. 

The second step involves the use of equation (1) which 
if we take logarithms and put a=F(/), becomes: 


log F—log 1.6=F (J) log 10t—-F(/) log J. (3) 
This equation can be simplified to: 
log log t’—F’ (7) (4) 


by ignoring the constant log 1.6, putting 10¢=¢’ and F'(/) 
log J=F’(J). Equation (4) in determinant form is: 


? 


—log 1 0 
log t’ 0 1 =0, (5) 
F’(J) 1 F(J) 


Finally the potential evapotranspiration (PE) is cal- 
culated according to: 


Ef=(PE) (6) 
which can be written: 
—m km log E 1 
n kn log f 1 =0 (7) 
0 kmn log (PE) m+n 


with m, and constants. 


3. CONSTRUCTING THE NOMOGRAM 


For the purpose of constructing the nomogram the de- 
terminant (5) must be multiplied by the non-zero de- 
terminant (see [6]) 
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a second nomogram is added to the first. Detern:inan 
| P2 Gs %T | #0 (8) (7) becomes with moduli 12.5 for the / line and 10 for the 
Ps Ys rs PE line: 

which gives: —2.5 12.5 log F 1 
—p, log E+ pez —q log E+ qe 10 50 log f 1 =0 (13) 

| 0 10 log PE 1 
+ pet (7) so the f line has a 50-em. cycle, while the distance between 
—r, log E+r the F and PE lines is 2.5 cm. and bet ween the PE andj 
r, log t’ +13 =0. (9) lines 10 cm. Since the # line now functions only as a 
(J) (J) pivot line it does not need to be marked. Markings for 


Inspection shows that it is permissible to put po=p~s= 
= qs=r"1=0, but pr, 72, and 7; must be non-zero. 
Therefore : 
—plogE 
pilogt’ 


in which 7, may be equal to 7;. By multiplying each row 
by the reciprocal of its element in the third column, ignor- 
ing subscripts, rearranging, and inserting the constant log 
1.6, we find the parametric equations: 


y¥:= (p/r) log 10¢ 

(p/r) (log log 1.6) 
F’ (J) 

TERS) 


=0 (10) 


(11) 


The nomogram can now be constructed and consists of 
two parallel logarithmic scales of equal modulus running 
in opposite direction and a third line not parallel to these. 
The slanted line is not straight because /'(/) is a poly- 
nomial. However, within the range of J’s for which 
Thornthwaite’s formula has proven to give adequate re- 
sults, the line is straight, except for slight curves below 
/J=10 and above J=160. Below /=5 the line approaches 
zero asymptotically toward (1/1.49239)(g/r) since 
by equation (2), /'(7)=0.49239 when J=0; and above 
J=160 the line curves similarly toward 7,=0. Between 
/J=10 and J=160 the slope of the J line is determined by 
t=26.5° C. on the ¢ scale and /’=13.5 em. on the £' scale. 
At this temperature # is 13.5 cm. regardless of index, 
For /'=1.6 em., equation (3) becomes: 


F(J) log 10t=F (J) log J (12) 


which, for /'(J)#0, gives log J=log 10¢. The ¢ scale 
can therefore be used to find the /’s on the sloping line by 
connecting 1.6 em. on the # scale with the different values 
(inem.) on the ¢ scale. 

For temperatures above 26.5° C. (79.7° F.) values of 
FE are given by Thornthwaite [4]. These are added to the 


F scale above 13.5 as indicated in figure 1. 

A convenient size of nomogram is obtained by making 
the distance between two vertical scales g/r=17 cm. and 
the length of one logarithmic cycle on the vertical lines 
p/r=12.5 em. 

In order to find the actual potential evapotranspiration 


°F. and inches can be added as desired to the ¢ and PE 
lines. 


4. ACCURACY AND PRACTICAL USE OF THE 
NOMOGRAM 


The nomogram gives results at least as accurate as those 
obtained by the old method and when used in routine pro- 
cedures has proven to be much faster, especially if one 
mean J value is established for any one station and a series 
of years of record is being processed. The deviations oe- 
curring in consequence of the use of a mean J value are 
negligible. 

Speed and accuracy both are increased even more if the 
nomogram is drawn on “hard board” or similar material, 
It can then be equipped with a system of pins and slots 
(for an illustration, see [7]). The J and the pivot line 
should be grooved. A pin is fastened at the appropriate 
place on the J line. A slotted transparent arm is made to 
slide over this pin and is connected via a sliding pin on 
the pivot line with another transparent arm. When the 
left arm is placed on the appropriate place on the ¢ line, 
the pin on the pivot line slides automatically into the right 
position and the right arm can be placed immediately on 
the proper place along the f line to obtain the PE. 

The f values depend on the month of the year and the 
latitude of the location. They must be obtained from 
tables as published by Thornthwaite ([4], p. 93) or 
Thornthwaite and Mather ([1], p. 98). Although it 
would be possible to arrange these tables in the form of 
graph, it would make a nomogram unnecessarily cluttered 
and awkward to handle. 

This nomogram, as any nomogram, is of principal ad- 
vantage if series of computations have to be made, sueh 
as series of years of monthly PE’s at one station or mea 
monthly PE’s at different stations. (In the latter case the 
stations should be grouped by latitude.) 

It is, therefore, convenient to mark the f values belong: 
ing to the latitude of the station, or stations, before taking 
a series of readings. As an example the mean temperatute 
for May in Bridgeton, N.J., is 17.5° C.; 7=58.3; /, taken 
from Thornthwaite’s table, is 1.23 at latitude 39° N.; and 
the PE is read to be 9.2 cm., as illustrated in figure 1. 


APPENDIX 


Table 1, as taken from Thornthwaite [4] presents the/ 
values for latitudes 50° N. to 50° S. for the 12 months of 
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A nomogram to determine the potential evapotranspiration according to Thornthwaite’s formula. 
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Tarte 1.—Mean possible duration of sunlight in the Northern and 
Southern Hemispheres expressed in units of 30 days of 12 hours each. 
(After Thornthwaite [4}.) 
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the year. For latitudes farther north or south the /’s of 
the 50th parallel should be used, presumably because the 
effect of longer daylight is offset by the diminishing effec. 
tiveness of radiation. 
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THE WEATHER AND CIRCULATION OF MARCH 1959 
Record Cold in Alaska But Mild Temperatures in the Remainder of the United States 


CARLOS R. DUNN 


Extended Forecast Section, U.S. Weather Bureau, Washingion, D.C. 


1. INTRODUCTION 


The departures from normal of circulation and tem- 
perature for March 1959 were rather small over most of 
the contiguous United States. Westerly flow, similar to 
the normal [6], prevailed in the mid-troposphere, and 
maritime polar air masses, with their moderating, spring- 
time effects, dominated the weather. However, the pre- 
cipitation was more abnormal. There were periodically 
heavy rain and snow falls as vigorous Lows developed in 
the Southwest and moved across eastern United States. 
Some regions received 200-300 percent of their normal 
precipitation. 

In sharp contrast to the mild weather in the contiguous 
States was the bitter cold throughout the State of Alaska. 
Long-period Alaskan records were broken, and the 
reported mean temperatures were more typical of mid- 
winter than early spring. For example, the monthly 
mean temperature at Barrow was a frigid 27° F. below 

Meanwhile, the Hawaiian Islands had unusually warm, 
dry weather. Honolulu, for example, reported a record 
high mean temperature and little precipitation. There 
were just 6 days with measurable rainfall, and the total 
for the month was only 7 percent of normal. 


2. THIRTY-DAY MEAN CIRCULATION 


The monthly mean circulation at 700 mb. for March 
1959 consisted of a large-amplitude, three-wave pattern 
in high latitudes and a smaller-amplitude, four-wave pat- 
tern in middle latitudes (fig. 1). These two wave trains 
were approximately 180° out of phase. Only the trough 
in Eurasia extended continuously across all latitudes, and 
even it had a marked slope from northeast to southwest. 
Associated with these truncated, out-of-phase troughs and 
ridges were several confluence zones with accompanying 
strong winds. The most prominent ones were in the At- 
lantic and eastern Pacific, where the height anomalies 
(fig. 1) indicate stronger-than-normal westerly geo- 
‘trophic winds. In fact, wind speed maxima were located 
in both of these regions, and mean speeds were more than 
$m.p.s. above normal (fig.2). The center in the Atlantic 
was & primary wind maximum, but the one in the eastern 
Pacific was essentially just an eastward extension of a 
larger primary wind maximum located in mid-ocean. All 
ofthese ind speed maxima were well east of their normal 


positions (indicated by arrowheads on heavy dashed lines 
in fig. 2). 

In Eurasia the mean winds were generally weak. Re- 
sidual blocking over Europe, which was associated with 
a split jet stream, produced a large area of subnormal 
wind over central Europe but supernormal wind to the 
north and south. One jet stream axis was depressed and 
flowed across northern Africa; the other was displaced 
north of the normal position, producing wind speeds up 
to 9 m.p.s. above normal between Scandinavia and Spitz- 
bergen. There were no organized jet streams at the 
700-mb. level over Asia. 

Returning to the 30-day mean circulation over the oceans 
(fig. 1), we note that the troughs as well as the primary 
wind maxima were displaced eastward from their normal 
positions. This eastward shift of the Asiatic coastal 
trough resulted in above normal heights and abnormally 
anticyclonic conditions over Japan and eastern Siberia. 
The cyclonic areas were also displaced eastward, produc- 
ing significant negative anomalies in the central Pacific 
Ocean and Alaska. 

The largest anomalies were in northern latitudes, where 
the 700-mb. heights were far enough above normal to re- 
sult in a High over Kamchatka, the approximate location 
of the normal Low, and also far enough below normal to 
produce a trough over Alaska, near the position of the 
normal ridge. This anomalous situation affected the 
downstream pattern, so that the trough normally found 
over eastern Canada was displaced approximately 25° of 
longitude eastward to southern Greenland, reducing 
heights in that area to 510 feet below normal, the largest 
anomaly in the Northern Hemisphere. 

In retrospect, it appears that the block over Europe and 
the unusually strong trough in western Siberia, through 
the flux of vorticity, molded the well-marked wave train 
downstream. In addition, the eastward displacement of 
waves was favored by stronger-than-norms| polar 
westerlies. 

In the United States 700-mb. height departures from 
normal were small (fig. 1). Over the western States posi- 
tive anomalies were associated with the flat ridge along 
the coast. In the eastern States there were weak negative 
anomalies over a broad band, typical of a flat mean trough. 
The latter feature suggests that cyclonic activity during 
March was uniformly distributed over the eastern two- 
thirds of the United States. 
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Figure 1.—Mean 760-mb. contours (solid) and height departures from normal (dotted), both in tens of feet, for March 1959. Mean 


ridge in the western and trough in the eastern United States were associated with stormy Weather east of the western mass. 


Abnormally deep trough brought cold weather to Alaska. 


3. CHANGE IN MEAN CIRCULATION 
FROM FEBRUARY TO MARCH 


The pronounced block, which was located over western 
Europe during February [2], was much weaker and far- 
ther east during March. This modification of the circula- 
tion produced over Great Britain the largest month-to- 
month height change in the Northern Hemisphere (—610 
ft.) and effected over Russia lesser height rises (fig. 3). 
This was associated with dry weather over Russia and 
_ much storminess over Great Britain. 

In other areas of the Northern Hemisphere the height 


changes (fig. 3) were similar to the 700-mb. height de 
partures from normal, which were referred to previous!) 
(fig. 1). Hence, in these areas the circulation anomalie 
for March can be looked upon largely as anomalous height 
changes from February to March. The rises over easter 
Siberia and the eastern Pacific, and the falls over Alask 
and central Pacific, all had their counterparts in the heigl! 
anomaly field. ‘To some extent the same was true in easter! 
North America, where heights rose in Canada and fel 
over the United States, resulting in weaker-than-norml 
west winds or height gradient from the southeaster 
United States to Hudson Bay. 
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Pisure 2--(A) Mean 700-mb. isotachs and (B) departures from 
monthly normal wind speeds, both in meters per second, for 
March 1959. Solid arrows in (A) indicate principal axes of 
maximum winds, and dashed arrows their normal March posi- 
tions. Regions with wind speeds greater than 15 m.p.s. and 
itomalies greater than 5 m.p.s. are stippled. Principal wind 


maxima were located over the oceans east of their normal 
sitions. 


4. RECURRENT CIRCULATION AND WEATHER 


The 5-day mean flow patterns during March resembled 
each other and were quite similar to the 30-day mean flow. 
This is illustrated by the geographical frequencies of 
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Ficure 3.—Difference between monthly mean 700-mb. height anom- 
aly for February and March 1959 (March minus February) in 
tens of feet. Largest falls were over Great Britain, Changes 
in many areas have striking resemblance to 700-mb. height 
anomalies for March (fig. 1). 


700-mb. 5-day mean troughs and ridges (fig. 4) and by 
their relationship to the 30-day mean chart (fig. 1). The 
bands of high frequency of troughs and ridges were nar- 
row, indicating a small variation in the longitudinal posi- 
tion of the troughs and ridges. One exception is the 
eastern United States, where the 5-day mean troughs were 
rather uniformly distributed from the Mississippi Valley 
tothe Atlantic Coast. 

It is particularly noteworthy that the positions of the 
5-day mean troughs and ridges clustered around their 
respective 30-day mean troughs and ridges. This means 
that the 30-day mean accurately portrayed the predomi- 
nant circulation regime of March and was not merely the 
average of several heterogeneous, short-lived circulations. 

The 30-day mean flow of small amplitude across North 
America was accompanied by zonal trajectories of the 
Highs for March (Chart IX in [4]). The continental 
polar anticyclones, which formed in western Canada, 
glanced east-southeastward and remained almost entirely 
in Canada. They affected only the northeastern United 
States where the monthly mean thickness (fig. 5) and 
surface temperatures (fig. 7) averaged below normal. The 
remainder of the United States was dominated by mari- 
time Pacific air masses, which also moved eastward across 
the country. A few anticyclonic centers had trajectories 
from the Pacific through the Northwest, but many of the 
Highs first appeared over the Plains and Central States 
and subsequently moved eastward into the Atlantic Ocean. 
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Fiecure 4.—Percent of the time that (A) troughs and (B) ridges 
on 5-day mean 700-mb. charts were located within 10° longitude 
intervals at latitudes from 20° N. to 70° N. for March 1959. The 
data were adjusted to an equivalent basis with 10° at 50° N. 
as the unit. Isoline interval is 20 percent. Areas with fre- 
quency greater than 20 percent are stippled; zero areas are 
hatched. Note high frequency of troughs and ridges near the 
locations of the 30-day mean troughs and ridges, respectively 


(fig. 1). 


These maritime air masses had a thickness or mean tem- 
perature of the layer from 700 to 1000 mb. that was above 
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Fieure 5.—Departure from monthly normal of the mean thick 
ness (700-1000 mb.) for March 1959 in tens of feet. Isoline 
interval is 50 feet. Below normal values are stippled. Patter 
over the United States bears a striking resemblance to the sur 
face temperature anomalies (fig. 7). Extreme negative de 
partures over Alaska were associated with record cold weather. 


normal for the northern areas of the United States. This 
was particularly true for Montana and the Dakotas, where 
in March continental polar air masses are required to pro- 
duce subnormal thickness. However, the same maritime 
Pacific air masses produced below normal] thickness in the 
southern States (fig. 5). 

The same was generally true for the accompanying sur- 
face temperatures (fig. 7), but this month an exception 
existed in the Far Northwest where subnormal surface 
temperatures were associated with above normal 30-day 
mean thicknesses. The coolest weather here occurred dur- 
ing the second and fourth weeks [7], when there was 
strong flow of abnormally cold air from the Arctic and 
Alaska (fig. 5) into the Northwest. These air masses had 
a short and rapid trajectory over the ocean, thus minimiz 
ing the diabatic heating, so that the air was still cool when 
it reached the northwestern United States. During the 
fourth week the cooling was further enhanced by cyclonit 
flow aloft. 

Along the California coast the weather was unusually 
warm. Both San Francisco and San Diego, which have 
long periods of record, reported new high mean temper 
tures for March. Furthermore, temperatures averaged 
above normal every day of the month at San Diego and 
were below normal only one day out of the month at St 
Francisco. Previously, several authors [1, 3] of thi 
series of articles have related the persistently above nor 
mal temperatures in recent years along the Califormi 
coast to abnormalities of the circulation and of sea-surfte 
temperature. This March the record high temperature 
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Fieure 6.—Tracks of selected major cyclones which produced intense storms over the United States during March. Circles are the 1200 
emt locations of cyclone on date indicated by number above position. Number below is central pressure of Low in tens and units of 
millibars. Bach week at least one severe storm moved out of the Southwest producing heavy snow and rain over large areas east of 


the Continental Divide. 


follow directly from the above normal heights at 700 mb., 
uid in addition, from the greater than normal north- 
easterly flow at 700 mb. and sea level, producing foehn 
winds for the area south of San Francisco (fig. 1, and 
Chart XT in [4]). 

The cyclones, with two exceptions, which had tracks 
over the United States also had origins there (Chart X 
in [4]). Several Lows formed in the Atlantic Seaboard 
States, but the storms which originated in the Southwest 
‘ontributed most to the weather over the United States. 
Each week at least one major cyclone formed over the 
southern Rocky Mountain States and produced severe 
weather along its eastward track through the central and 
‘astern United States. A detailed descripition of the 
tflects of these storms is given in [7]. 

During the first week of March a major storm marked 
by heavy rains, damaging winds, blowing dust, and drift- 
ig snow moved from Colorado northeastward across the 
“reat Lakes into Canada (fig. 6A). The weather during 
thesecond week was influenced mainly by two disturbances 


which reached strong storm intensity as they moved north- 
eastward (fig. 6B). The first storm, 11th to 13th, moved 
erratically from the lower Mississippi Valley to a posi- 
tion just off the north Atlantic coast and finally turned 
northward across eastern Maine. This storm produced 
heavy rains in the South and heavy snows in the North- 
east, as high winds raked the entire east coast. More 
heavy rains in the South and a band of heavy drifting 
snow from Nebraska to Michigan resulted from the second 
storm of that week (14th to 16th) which moved from the 
lower Great Plains across the Great Lakes. During the 
third week the incidence of Southwest storms continued, 
but the cyclone which formed over Utah on the 19th moved 
in an east-southeasterly direction with its path traversing 
the Gulf Coast States (fig. 6C). In addition, a cyclone 
formed over Lake Ontario and slowly traced an erratic 
track across the Maritime Provinces of Canada (fig. 6C), 
creating stormy conditions in New England on March 21. 
The mild, sunny weather east of the Rocky Mountains 
during the first part of the fourth week of March was 
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IF'icure 7.—Departure of average surface temperature from normal 
(° F.) for March 1959. Largest anomalies were in the northern 
Plains and Rocky Mountain States, but the positive anomalies 
along the California coast were associated with record tem- 


peratures (from [7]). MARCH 1959 


Figure 10.—Number of days in March 1959 with fronts of an 
type within unit squares (with sides approximately 500 miles) 
All frontal positions are taken from Daily Weather Map, 18 
p.m. EST. Areas with 15 or more days with fronts are stippled 
Active fronts were frequently located in Florida, Central Plains 
and the Northwest. 


interrupted by stormy weather as still another Texas dis 
turbance moved northeastward on the 25th to 27th frou 
the Southwest across New Jersey (fig. 6D). 

The severity of the weather associated with thes 
“Southwest Lows” is attested by the summary of tle 
monthly weather at Dubuque, Iowa, which felt the ful 
brunt of several of these storms. It was reported thi’ 
the precipitation for March was the heaviest since tle 
beginning of records in 1851. The snowstorm of the 4) 
and 5th brought the greatest amount of snow ever It 
corded in one 24-hour period for any month. The 6.) 
inches of moisture for the month exceeds by 1.48 inehe 
the previous wettest March recorded in 1852. The total 
monthly snowfall of 30.2 inches is the greatest for ai) 
March and the greatest for any month of the year sine 
January 1929. Similar reports of record-breaking weath! 
came from stations in Nebraska and Wisconsin. 

The bulk of the precipitation over the United State 
during March (figs. 8, 9) was produced by the majet 
storms whose tracks are shown in figure 6. Howevel 
additional rainfall did occur in the Southeast, associate! 
with fronts which were prevalent over Florida and tlt 
Gulf of Mexico (fig. 10). Some of the precipitation 
cold-frontal in nature, but the majority was produced }) 
incipient cyclone waves and over-running of the col 
air by tropical air masses. Tampa reported the welte® 
for March on record; other Florida stations broke yaridl 

€ e so eastern, nor ’ 

but dry in the north- precipitation records. 

In the Northwest above normal precipitation occw™ 


ern Plains, northern Rocky Mountains, and the Southwest i 
(from [7]). with faster-than-normal west-southwesterly flow at 


Ficure 8.—Total precipitation (inches) for March 1959. Record 
amounts fell in Florida (from [7]). 
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the mountain ranges (fig. 1), which intensified orographic 
lifting. Also favorable were the frequent invasions of 
maritime Pacific fronts (fig. 10) and cyclonic flow aloft. 
(See Duily Weather Map [5].) 

That precipitation which fell east of the Rocky Moun- 
tains occurred under cyclonically curved 700-mb. mean 
fow and weak easterly anomalous flow in the Kansas- 
Missouri-Illinois area; both types of flow have been re- 
lated empirically to vertical motion and precipitation. 
However, the explanation for the precipitation which oc- 
curred in the easterly flow north of the daily Lows, along 
the eastern slopes of the Rocky Mountains in Colorado 
and Wyoming, is not easily extracted from the monthly 
mean circulation at the 700-mb. level. The vertical mo- 
tion and precipitation in these central Rocky Mountain 
States were enhanced by orographic effects not reflected 
in the upper-level circulation. Also, the precipitation 
occurred during a few short periods which were not typi- 
cal of the predominant pattern of the region and conse- 
quently would not be reflected in the 30-day mean 
circulation. 


5. RELATION OF PRECIPITATION TO VERTICAL 
MOTION 


Of course, subjective evaluation of the monthly mean 
flow in terms of the average vertical motion or precipita- 
tion is rather crude, but until recently, it has not been 
practical to objectively estimate the vertical motion. Now, 
the availability of operational baroclinic models and elec- 
tronic computers enables us to compute, objectively, the 
vertical motion for a midtropospheric level. The National 
Meteorological Center of the U.S. Weather Bureau at 
Suitland, Md., using a two-level baroclinic model, esti- 
mates twice daily the concurrent vertical motion ‘at 600 
mb. from initial (observed) data. 

The 30-day mean of these values for March 1959 (fig. 
ll) is a logical pattern roughly consistent with the con- 
comitant mean 700-mb. circulation. In general in figure 
ll, descending motion (negative values) is found west of 
the troughs in the northwesterly flow and ascending mo- 
tion (positive values) east of the troughs in the southerly 
flow. Exceptions are the negative values just off the north- 
eastern coast of the United States, under and east of the 
mean trough (fig. 1,11). This, at least at first thought, 
does not fit the accepted relationship between horizontal 
circulation pattern and associated vertical motion, but 
further inspection of figure 1 reveals that subnormal 700- 
mb. heights and marked anomalous, northerly flow ex- 
tended from eastern Canada to New England, suggesting 
that cold subsiding air masses frequently invaded the area 
in question, off the northeastern coast of the United States. 

Portions of the 30-day mean vertical motion were cor- 
related with the total monthly precipitation (fig. 8). The 
descending motion over the extreme Southwest and the 
uorthern Plains States was related to the driest areas this 

Prescott, Ariz., which had no precipitation, and 
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FieurE 11.—30-day mean vertical motion at the 600-mb. level in 
millimeters per second for March 1959. Map was obtained by 
averaging the 60 available twice-daily values computed from 
baroclinic model and observed data. Isoline interval is 2 milli- 
meters per second. Absolute values greater than 2 are stippled. 
Descending motion (negative) and ascending motion (positive) 
tend to be located east of the 30-day ridges and troughs, 
respectively. 


Ficure 12.—30-day mean of only the ascending (positive) daily 
vertical motion for March referred to in legend of figure 11. 
Isoline interval is 2 millimeters per second. Values greater 
than 4 are stippled. Pattern resembles the total monthly pre- 
cipitation (fig. 8). 


Los Angeles, Calif., reported record dryness for March. 
Helena, Mont., had the driest March since 1881. The 
ascending motion over the Far Northwest and the States 
bordering the Gulf of Mexico fit the precipitation pattern 
quite well. However, pronounced discrepancies existed in 
the Nebraska-Kansas-Iowa and Northeast areas where 
sizable amounts of precipitation occurred in regions of 
little or no mean ascending motion. This could result 
from the averaging process where large ascending motions 
during short periods of storminess would be counteracted 
by long periods of subsidence. Therefore, the mean of 
only the positive values (ascending motion) was obtained 
(fig. 12). As expected, there was considerable ascending 
motion in these wet regions. Apparently, west of the 
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mean troughs, on the western fringe of the precipitation 
areas, there were periods of both ascending and descending 
motion associated with the lifting and precipitation north 
and east of the daily storms and the subsidence in the 
cool air behind these same storms as they moved eastward. 
The availability of moisture, which is an additional re- 
quirement for precipitation, has not been considered in this 
brief investigation, and should be included in any further 
study. 


6. ALASKAN COLD AND HAWAIIAN WARMTH 


Most of Alaska experienced a record-breaking cold 
March (table 1). Barrow, Fairbanks, and Nome, well 
distributed stations with long periods of record, had the 
coldest March on record and Anchorage the second coldest. 
Departures from normal of the monthly mean temper- 
atures were as large as two or three standard deviations, 
a rather rare event. Fairbanks, with 16° F. below normal, 
had the most extreme anomaly, but Barrow, with a frigid 
mean temperature of 27° F. below zero, suffered the coldest 
weather in the absolute sense. At the latter city the maxi- 
mum daily temperature did not exceed —5° F., and the 
minimum was always colder than —21° F. 

This cold over Alaska was produced by stronger than 
normal flow between the mean ridge over eastern Siberia 
and the trough over Alaska (fig. 1), which continually 
advected cold air masses from the Arctic Basin into 
Alaska. This was a very persistent pattern, as illustrated 
by the 5-day mean trough and ridge frequencies (fig. 4) 
and produced over Alaska the coldest mean thickness 
anomaly for the Northern Hemisphere (fig. 5). 

While the Alaskans were shivering in record-breaking 
cold weather, the Hawaiians basked in record-smashing 
warmth (table 1). At Honolulu, the daily maxima aver- 
aged 79° F. and the minima 71° F. for the month of March. 
Of more significance was the monthly mean temperature 
of 74.7° F., the absolute highest since records began in 
1905. This month’s mean temperature represents a de- 
parture of +2.5° F., which is large for this subtropical 
station that has a standard deviation of only 1.0° F. 
Similar weather was experienced at Lihue, Kauai, where 
a record monthly mean temperature of 74.3° F. produced 
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an anomaly of +3.7° F., almost 3 standard deviation; 
above normal. 

This warmth was associated with a deficit of precipita. 
tion at both stations. Honolulu reported only 0.17 ineh, 
2.13 inches below normal, and Lihue had 1.37 inches, 27 
inches less than normal. 


1,—Alaskan and Hawaiian surface temperatures (° PF.) fy 


March 1959 
Monthly Standard! Yegr 
Station mean | Normal | Anomaly |deviation | records 
Start 
Alaska 
14.0 24.8 —10.8 4.6 1916 
*—6.6 9.0 —15.6 6.0 1906 
33.8 32.8 +1.0 2.8!) 1915 
Hawaiian Islands 
*74.7 72.2 +2.5 1.0 | 1905 
"74.3 70.6 +3.7 13] 1905 
*Record for March. 


The dry and extremely warm weather occurred with 
a stronger than normal subtropical ridge in the eastem 
and central Pacific. Over the Hawaiian Islands 700-mb, 
mean heights (fig. 1), 1,000-700-mb. mean thicknesses, 
and sea level mean pressures were all above normal dur 
ing March. 
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CORRESPONDENCE 
(Continued from p. 106) 


We were not aware of a similar modulated technique in 
use by the Weather Bureau. 

We are aware that saturation thresholds are somewhat 
difficult to work with, but for conditions of level plot land, 
we assume that a saturation point exists and that any 
precipitation in excess of this value is runoff and hence 
lost as far as crop use is concerned. 

During the past year we have had the opportunity to 
compare atmometers with the Class A pan, and the buried 


4-foot pan. These observations have not caused us! 
alter our point of view on the disadvantage of pans, ge 
erally. Frost is a serious limitation to the use of atmol 
eters. (In this connection we have been experimentilf 
with an instrument which withstands 7 degrees of frost. 
Our co-workers in Canada have noticed time trends» 
Bellani plate atmometers only when air enters the cup 
when the surface color or porosity changes because d 
dirt or accumulated salts. 
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